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Abstract 

One of the key problems related to the Bak-Sneppen evolution 
model on the circle is to compute the limit distribution of the fitness 
at a fixed observation vertex in the stationary regime, as the size of 
the system tends to infinity. Simulations in [I] and suggest that 
this limit distribution is uniform on (/, 1), for some / ~ 2/3. In this 
paper we prove that the mean of the fitness in the stationary regime 
is bounded away from 1, uniformly in the size of the system, thereby 
establishing the non-triviality of the limit behaviour. 

The Bak-Sneppen dynamics can easily be defined on any finite con- 
nected graph. We also present a generalisation of the phase-transition 
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result in the context of an increasing sequence of such graphs. This 
generalisation covers the multi-dimentional Bak-Sneppen model as 
well as the Bak-Sneppen model on a tree. 

Our proofs are based on a 'self-similar' graphical representation of 
the avalanches. 

1 Introduction 

The Bak-Sneppen model, introduced in [2], has received a lot of attention 
in the literature, see for instance [I], [I] and [Zj. In [I], it is described how 
Bak and Sneppen were looking for a simple model which was supposed to 
exhibit evolutionary behaviour, and which was also supposed to fall into the 
class of processes showing self-organised critical behaviour. For physicists, 
self-organised critical behaviour refers to power law decay of temporal and 
spatial quantities. After a number of attempts, Bak and Sneppen arrived at 
the following simple process. 

Consider a system with N species. These species are represented by iV 
vertices on a circle, evenly spaced, say. Now each of these species is assigned a 
so called fitness, a number between and 1. The higher the fitness, the better 
chance of surviving the species has. The dynamics of evolution is modelled 
as follows. Every discrete time step, we choose the vertex with minimal 
fitness, and we think of the corresponding species as disappearing completely. 
This species is then replaced by a new one, with a fresh and independent 
fitness, uniformly distributed on [0,1]. So far, the dynamics does not have 
any interaction between the species, and does not result in an interesting 
process. Indeed, if we only replace the species with the lowest fitness, then 
it is easy to see that the system converges to a situation with all fitnesses 
equal to 1. Interaction is introduced by also replacing the two neighbours 
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of the vertex with lowest fitness by new species with independent fitnesses. 
This interaction represents co-evolution of related species. This neighbour 
interaction makes the model also very interesting from a mathematical point 
of view. 

It is simple to run this model on a computer. Simulations then suggest 
the following behaviour, for large N (see |1] and pQ for simulation results). 
It appears that the one-dimensional marginals are uniform (in the limit for 
iV — > oo) on (/, 1) for some / whose numerical value is supposed to be close to 
2/3. This threshold value / is the basis for self-organised critical behaviour, 
according to [2], [I] and jl], as follows. Since in the limit there is no mass 
below /, one can look at so called avalanches of fitnesses below this threshold: 
starting the counting at the moment all fitnesses are above / and finishing 
the counting at the first next moment all fitnesses are above / again. The 
random number of updates, for instance, counted this way, is supposed to 
follow a power law. For this to make sense, the conjecture had better be 
true. This we have not been able to prove. However, we have been able 
to prove the weaker result that the mean average fitness in the stationary 
regime is bounded away from 1, uniformly in the number of vertices. Hence, 
in the limit, there is probability mass of the fitnesses below 1. Since it is not 
difficult to show that in the limit there can be no probability mass of the 
fitnesses below 1/3, our results establish that the one- dimensional marginals 
do not become trivial as iV — > oo. For a similar result in a discrete version 
of this process, see [6J. 

Let Fn be the distribution function of the one-dimensional marginal in 
the stationary regime, in the system with N vertices. We will prove the 
following result. 

Theorem 1.1 If q < 1 is close enough to 1, then there exists c q > 0, inde- 
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pendent of N, such that 

F N (q) > c q . (1.1) 

In Section 0] and 5, we define the Bak-Sneppen dynamics on arbitrary 
finite connected graph, and present a generalisation of our results for an 
increasing sequences of such graphs. 

We remark that corresponding results for a mean-field version of the 
model are quite simple to obtain, see e.g. jHj and In the mean-field 
case, it is possible to prove, using essentially only combinatorics, that the 
one-dimensional marginals do indeed converge to a uniform distribution on 
(c, 1), for some constant c which depends on the characteristics of the model. 

In the next section we prepare for the proof of Theorem II. 1[ by introduc- 
ing the notion of an avalanche, and establishing some monotonicity properties 
of the avalanches. The proof itself can be found in Section 3. 

2 The self-similar graphical representation 

Let A(iV) = {— N + 1, . . . , —1, 0} index the set of N vertices on the circle, 
so that and —N + 1 are neighbours. We use negative indices to simplify 
notation in the future. We say that in the time interval [n, n+d], an avalanche 
from threshold q G [0, 1] (also referred to as a g-avalanche) with origin at 
x G A (AT) and duration d > 1 occurs if at time n, x is the vertex with 
minimal fitness, above threshold q, and n + d is the first moment after n with 
all fitnesses again above q. The range set of the g-avalanche is the collection 
of vertices updated during the avalanche, and the range of the g-avalanche 
is the number of different vertices in the range set. Note that, according to 
this definition, if at times n and n + 1 all the fitnesses are above q, then in 
the time interval [n, n + 1] an avalanche of range 3 and duration 1 occurs, 
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even though there were no fitnesses below q. 

The fitnesses of the vertices are random variables with values in [0, 1] 
and we update them according to the uniform distribution on [0, 1]. For 
computational reasons however, it is convenient for the fitnesses to have 
values in [0, oo] and to update them according to the exponential distribution 
with parameter 1, say. In this new setup a threshold b corresponds to the 
threshold q = 1 — e~ b in the old setup. 

Suppose that a 6-avalanche starts at time 0, with the origin at the vertex 
0, so that the vertex and its two neighbours are updated. We can now 
graphically illustrate the 6-avalanche on A(iV) x M + (space x fitness) as 
follows. Look for the vertex with minimal fitness, and call this vertex x. 
(Note that x must be the vertex or one of its two neighbours.) Suppose 
that the fitness of x is equal to s < b. We then continue updating according 
to the appropriate rules, and wait until all fitnesses are above the threshold 
s. This in itself constitutes an s-avalanche, starting at x. We denote by 
Cn{x, s) the set of vertices involved in this s-avalanche. In the graphical 
representation, we draw an arrow from the space-fitness point (x, s) to the 
space-fitness points (i, s), for all i G £jv(x, s )- 

After the s-avalanche has ended, the new fitnesses of all vertices involved 
in this avalanche are i.i.d. and exponentially distributed on [s, oo), due to 
the lack of memory property of the exponential distribution. We can now 
look for the minimal fitness among all vertices in £n(x,s). If this minimal 
fitness is above b, then the 6-avalanche has stopped. If this minimal fitness 
is equal to t, where s < t < b, and is associated with the vertex y, say, then 
we start, as before, a t-avalanche with origin y. We continue updating until 
all fitnesses are above t. If £n(v, t) denotes the set of vertices involved in this 
t-avalanche, then we draw an arrow in the graphical representation from the 
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space-fitness point (y,t) to all space-fitness points (i,t), for % G £/v(?/, t). We 
continue in the obvious way. Under the assumption that all avalanches are 
finite, this process will stop a.s. as soon as all fitnesses are above b. The 
idea of avalanches which form a hierarchical structure of subavalanches is 
also mentioned in in a slightly different context. 

This graphical representation, denoted by GRn, is a random graph on 
the space-fitness diagram A (AT) x R + . We can describe it more formally as 
follows. 

Let {n fc } fcg A(iv) be a collection of independent homogeneous Poisson pro- 
cesses. For each process IT^. we perform the following procedure. At the j-th 
arrival r kjj of U k , we draw a pair (£ N (k,T kJ ),r) N (T k j)J , where £ N (k, r kij ) 
is distributed as the range set, f)N(T k ,j) as the duration of a typical r k ,j- 
avalanche, with origin at k. We draw arrows in A(N) x M + from (k,r k j) to 
{y, r kt j), for all y G £at(£;, T k j). For any t\ < t 2 we say that (x, ti) is connected 
to (x, t 2 ) by a time segment. A path is a sequence (xo, s ), . . . , (x n , s n ) of 
points in A(N) x M + such that every pair (xj,Sj), (xj+i, s_j+i) is connected 
by either by a time segment or an arrow. For any x, y G Z, and t\ <t 2 a M, 
write (x, ti)-^(?/, t 2 ) in GRn, if there exists a path from (x, ti) to (y,t 2 ). See 
Figure n for an illustration. 

For any b > the range set £at(0,6) of a 6-avalanche with origin at 
consists of all vertices x such that (0,0)-w(x, b) in GRn, and the duration, 
denoted by rj N {b), of this avalanche is the sum of fj N {r X j) over all t x j < b, 
such that (0, 0)-^(x, t x j) in GRn- 

The graphical representation provides us with the following monotonicity 
properties. For any A C A(N) and t,s > we denote by £jy (s) the 
random set of vertices s G A (A/ - ) such that there exists y = y(x) G A, and 
(y, t)^(x, t+s) in GRn- Similarly, for any A C A (AT) and t, s > 0, we denote 
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Figure 1: The graphical representation GRn, where, for instance, (fc, — 2,t + s) 

in GRn- 

by r)x (s) the sum of t) n {t x ^) over all r x j < t + s, such that there exists 
y = y(x) G A, and (y,t)~^(x,r x> i) in GRn- Then for any A C 5 C A(iV), 
< si < s 2 , and t > 

In particular, for ^(0, 6) = 4 { ~ 1,0 '~ W+1} ' 0) (&) and 77^(6) = r/"- 1 ' '"^ 11 ' ^), 
we have 

£at(0A) ceAr(0,6 2 ), ^iv(fci) < tjjvfo), if6i<& 2 - ( 2 - 2 ) 

The last inequality allows us to couple two copies GRn and Gi?^ of the 
graphical representation in such a way that for any k G A(iV), IL^-) restricted 
to [0, 6/2] is the same as n' fc (-) restricted to [6/2, 6], and fj{r k j) < f)'{r' k .), for 
r fc j - G [0,6/2] and 7j£ • G [6/2,6]. See Figure El for an illustration of the cou- 
pled G-Rat and GR' N . Note that we do not claim that the two copies together 
yield a realisation of the evolution of the process from to 6. Nevertheless, 
this coupling gives us, for any 6 > 

V ^°\b/2)<4 MN)m (b/2), 
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Figure 2: The lower part of GRn and the upper part of GR' N coupled together, 
and hence 



2E( V ^ NW2) < 



\b/2))> E {rj^°\b)). (2.3) 

3 Proof of Theorem 11.11 

An important step in the proof of Theorem 11.11 is the following lemma, esti- 
mating the probability that an avalanche has range N, uniformly in N. For 
any b G define P/v(6) as the probability that an avalanche has range N. 

Lemma 3.1 If b is large enough, then P/v(5) > 1/2, uniformly in N. 
PROOF: For any x G A C A(N), denote by £(x, A) the left corner of A with 



respect to x, 



(x, A) = min jfc G (— oo, x] [k, x] C Amod iv|. 



and write In(s) := £ (0, £n(0, s)) for the leftmost vertex involved in an s- 
avalanche with the origin at 0. We will have proved the lemma if we show 
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that In is explosive in the sense that there exists a < < oo such that 
for any < i < i max = [log 3/2 (iV - 1)] G N, 

P (l N (boo) < max { - (I)* , -N + l}) > \ + (§) i+1 , (3.1) 

where the 'max' actually only works at % = i m ax- Indeed, (13. 1|) implies that 
-P/v(&oo) > |. In order to achieve this, choose a constant bo > 17. Define a 
converging sequence of thresholds bi, b 2 , b 3 , . . . as follows: 

h= k-l+(f)*&0, Z>1, 
fooo = lim bi = 4b Q . 

i— >oo 

Observe that due to the monotonicity property (j2.1|) . it suffices to prove that 
for all i G [0, i max ] 

P (l N {h) < max { - (§)' , -N + l}) > \ + ( 1 i y+ 1 . (3.2) 
We proceed by induction. First note that 

P MM <-!) = !■ ( 3 - 3 ) 
Next, suppose that (|3.2j) holds for some i G [0, i max — 1]. Observe that 

l N ibi)<-%\ 
< 3x G [— (f) 4 , -| (I) 4 ] n A(iV), 3^- G U x n [6,, 6 m ), > 
such that £ ^a;, ^(x, r fcj ) j < x — (|)* 

implies 

{ < max {- (|) m , -N + l} } . 

See Figure|3]for an illustration. Hence to finish the inductive step, it suffices 
to show that for all i G [0, i max — 1], 

( v.x g [- (§)* , -| (I) 4 ] n A(iV), Vr XJ g n x n [6,, fe m ), 

y l(x,£ N {x,T k j)} > X- (§)' 

(3.4) 
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Figure 3: Illustration of the induction step. 

Since the £ (x, iu(x, r^j) Vs are independent and since (due to the mono- 
tonicity property) for any x G A (AT) and t x j > bi, 

P (£(x,i N (x,T xJ )) >x- < P (l N {bi) > - (I)') < |, 

the points 

*e[-(f)\-§(f)*]nAW 
constitute a thinning of the Poisson process (J IT^ with 

deleting probability at most 1/2. Thus the points in ()3.5|) contain a Poisson 
process of intensity at least 

\ ■ ( intensity of \J n x J 

Observe that the event in (|3.4j) implies that the process (|3.5|) has no arrivals 
between time 6, and bi + i, a time interval of length &o(|) 4 ; and hence it has 
probability at most 

-P{-| (§)' bo (I)*} = exp{-| (§)<} < (iy +2 , i G N, 
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since b > 17. So we have ()3.4|) and the proof is complete. 



□ 



Proof of Theorem ll.lt Suppose we start from an i.i.d. uniform distri- 
bution above the threshold q < 1, which we assume however is so close to 
one that in the model where we update fitnesses according to the exponen- 
tial distribution, we would have Pjj(b/2) > 1/2, uniformly in N. (Recall 
that b and q are related via q — 1 — e~ b .) We will define the dynamics via 
the following independent sequences of i.i.d. random variables. Fix some 
q' E (q, 1). Let, for i = 1,2,3, W = (Uj)j e j$ be a sequence of i.i.d. random 
variables uniformly distributed on [q',1]. We will use W to construct the 
dynamics above threshold q'. Let, for i = 1,2, 3, V 1 = (V?)j 6 n be a sequence 
of i.i.d. random variables uniformly distributed on [0,g']. We will use V 1 to 
construct the dynamics below threshold q'. Let, for i = 1,2,3, S l = (Sj)j G n 
be a sequence of i.i.d. Bernoulli distributed random variables taking the value 
1 with probability l — q' and the value with probability q'. We will use S i to 
choose between W and V\ and store the result in the sequence Q % = (G*)j £ N, 



It is clear that for i = 1,2, 3, the sequence Q % consists of i.i.d. random variables 
uniformly distributed on [0, 1]. We will use the beginning of the sequence Q 2 
(this choice is arbitrary) to assign the initial fitnesses to the vertices. At time 
n E [0, N — 1], we use the random variable G 2 n to assign an initial fitness to 
the vertex with number — n. Now we are ready to define the dynamics. At 
every time n > N we choose the vertex k, say, with minimal fitness at time 
n — 1, and we assign to (k — 1, k, k + 1)( mod N) the triple (G^, G\, G^). 



i = 



1,2,3: 
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Let fo(n) denote the fitness of a fixed observation vertex at time n. Let 
j n be the moment that this vertex received its current fitness, i.e., 

j n = min \j<n\ f (j) = / (n)|. 

Let i(n) be the number of the sequence providing this value, i.e. 

/o(n) = Gf \ 

Observe that (since jj n = j n ) i(n) = i(j n ). We say that n is q' -good, if during 

[j n , n] the minimal fitness of all vertices is always less than q'. Then, for any 

time fiGN and any q" G (q', 1), we have 

n-l 

P(fo(n) < q") = 5>(/o(n) < q", Jn = J 

3=0 
3 n-l 



i=l j=0 
3 n-l 

> XH^ P ( G 1 < = = i, n is g' - good 

8=1 j = 

3 n-l 

> Y Y P { U d < q "> i n = 3> W = h n is q' - g° od ) • 



i=\ j=0 



(3.6) 

Now observe that for any j G [0, n], the event |j n = j, = i, n is q' — 
good f is measurable with respect to 



a{ll}„ Vj„ S)„ Sp,, Vp„ < f < j, j < j" < n, i = 1, 2, 3}. 

Hence, for any j G [0, n], the events |Z7j < g"j and |j„ = j, = 
i, n is g' — goodj are independent, and we can continue the estimate (|3.6|) 
as 

3 n-l 

= Y Y P { U i < ^") p (^ = = *, n is g' - good 



i=l i=0 

3 n— 1 

'rEf YY P n = = *, n is g 7 - goodj 

i=l j=0 

g"-<?' 



j3^P(^n is g' — goodj . 
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It remains to estimate the probability that n is g'-good from below, uniformly 
in n and N. 

Define a sequence (r^)j e ^ of stopping times, with respect to the natural 
filtration, as follows: r = 0, and for any j G N, Tj +1 is the end of the first 
g-avalanche of range iV after Tj. For any j G N, we call the time interval 
I?* = [Tj,Tj+i) the j-th period. It is clear that at every t^ the fitnesses 
are i.i.d. and uniformly distributed above the threshold q. Thus the period 
lengths are i.i.d. random variables. For any time n, we denote by Tji n ) the 
maximal Tj such that Tj < n, i.e., n G [Tj( n ), T j(n)+i), and we say that n is 
q'-nice, if during [rj(„),n] the minimal fitness is always less than q' . For any 
n G N, if n is g'-nice then n is g'-good. Indeed, suppose that n is g'-nice. 
If jn > T j(n), then n is clearly g'-good. Suppose j n < Tj( n y Since during 
the g-avalanche of range iV of the previous period, every vertex has been 
updated, j n belongs to this g-avalanche, and hence the minimal fitness at the 
time interval \j n ,Tj(n)) is always less than g < g', and n is g'-good. Thus it 
suffices to show that the probability that time n is g'-nice, is bounded away 
from zero, uniformly in iV and n. 

A period can be decomposed into two parts: the duration of the avalanche 
of range N, and the waiting time until this avalanche. We denote by a 
typical waiting time before the avalanche of range N, and by An the duration 
of this avalanche. During a g-avalanche, the minimal fitness is always at most 
g < q' . Hence, if in the i-th period, the waiting time satisfies = 0, 
and in addition there is at least one vertex at time with fitness between 
q and q', then any time n within the i-th period is g'-nice. The event that 
there is such a vertex with fitness between g and g' is independent of Wn 
and A N associated to that period, and has probability 




uniformly in N. 
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Hence, in the stationary regime with N vertices, we can write (using alter- 
nating renewal process theory), 

P{n is (/-nice) > p\P{n is in a period with Wn = 0) 

for n — > oo. At this point, we switch from q to b, since we use results from 
the previous section. 

Since P(W N = 0) = P N (b) > P N {b/2) > 1/2, uniformly in N, it suffices 
to prove that there exists a constant < c(b) < oo, independent of N, such 
that 

E(W N ) < c{b)E{A N ). (3.8) 
Denote by Y N the number of 6-avalanches preceding the 6-avalanche of range 
N. Every avalanche has range N with probability Pat(6), independently of all 
other avalanches. Hence Yjv+1 has a geometrical distribution with parameter 
Pn (b), and we have 

Let (Z^)^ be an i.i.d. sequence of random variables distributed as the 
duration of a typical 6-avalanche, conditioned on its range being smaller 
than N. Then we can use Y N of those avalanches to obtain Wn, i.e. 

w N = z? + -.- + z? N . 

In words, at the beginning of a new avalanche, we first decide (with the 
correct probability) whether or not the avalanche has range N. If not, we 
choose one, conditioned on its range being smaller than N, and the resulting 
duration is the next Zf . Since due to the construction, Y N is independent 
of the sequence (Z^) i£N , we have that 

E(W N ) =E(Z» + ... + Z» N ) 

= E{Y N )E{Z?) = (p^y - l) E{Z»). 
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We will now estimate E(Z±) from above and E(An) from below (recall that 
the pair {£,N{b),r} N (b)) represents the range set and duration of a 6-avalanche) 

E(Z?) =E(r) N (b)\\£ N (b)\<N-l) 

^ kP( VN (b) =k,\£ N (b)\ <N-1) 
U P(Mb)\<N-l) 

y> kP( VN (b) = k) 

~ P(\Ub)\ < N - 1) 
<E(ri N (b))/p(\Z N (b)\<N-l 

E(A N ) =E(ri N {b)\\£ N (b)\ = N 



>E(p N (b)l{\£ N (b)\ = N}^ 
>E(r] N (b)l{\£ N (b/2)\=N}^ 
> E[ V ^ NW2 \b/2)l{\Ub/2)\ = N}^ 

and, since rf^ N ^' b ^ (b/2) and 1{|£jv(&/2)| = N} are independent, this is 
equal to 

^(^ A(JV) ' V2) (6/2))P(|^(6/2)| =N)> ^E(rj N (b))p N (b/2), 

where the last inequality follows from (|2.3J1 . Combining the estimates of 
E(Z^) and E(A N ), we have 

Since P/v(&) > Pat (6/2) > 1/2, the above inequality together with (jH.fjjl gives 
us 

£(Wat) < (p^y - l) rzj^j • jv§72)£(^v) 
Thus we have (|H.8jl and the theorem. 

□ 
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4 Extension to general graphs 

Let G be a finite connected graph. One can define on G a Bak-Sneppen 
process in the following way. We call two vertices neighbours if they are 
connected by a bond in G. Every vertex of G accomodates a random variable 
(the fitness), with value in [0, 1]. At the initial moment, all the fitnesses are 
i.i.d. and uniformly distributed on [0, 1]. Every discrete time step we choose 
a vertex with minimal fitness, and replace it, together with the fitnesses of all 
its neighbours, by new independent fitnesses, uniformly distributed in [0, 1]. 
We will give bounds for the mean of the fitness in the stationary regime. 
These bounds will depend on the local geometrical structure of G, but will 
be independent of the number of vertices in G. 

In order to state our main result, here follows some notation. Let V G 
denote the set of vertices of G. For any two vertices x, y G V G) we denote by 
p G (x, y) the distance between them, i.e. the number of bonds in the shortest 
path between x and y. Then for any vertex x G V G and k G N, we can define 
the ball 

B%(k) = {yeV G \p(y,x)<k}, 

and the sphere 

S G (k) = {yE V G \p G (y,x) = k}. 

Observe that since G is finite, we have, for sufficiently large k = k{G), that 
B G (k) = V G and S G (k) = 0, for any x G V G . 

For any k G N denote by m G (k) the number of vertices in the smallest 
ball of radius k and denote by M G (k) the number of vertices in the largest 
sphere of radius k, i.e. 

iTi G (k) = min \B G (k)\, 

x£Vg 

M G (k) =max\S G (k)\. 

x€Vg 
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Consider an increasing sequence of radii 

% E N. (4.1) 

For any vertex x of G we denote by Fq the distribution function of the 
fitness at x in the stationary regime. In the following two theorems we will 
establish an analogue of the phase-transition result for an infinite collection 
of finite connected graphs. 

Theorem 4.1 Let Q be an infinite collection of finite connected graphs such 
that for some 60 < b < oo, 

|(|) m G (h/3j) > log (M G ([2rV3l)), uniformly in G E Q,i E N. 

(4.2) 

Then there exists q, > and % E (0, 1) such that 

^g(?&) > Cb, uniformly in G E Q, x E G. 

Theorem 4.2 Let Q be an infinite collection of finite connected graphs of 
uniformly bounded degree, i.e. there exists a constant K E N such that 

max |-B G (1)| < K, uniformly in G E Q. (4.3) 

xEVq 

Then for any sequence (G n ) ng N C Q such that 

|V G J — > oo, as n — > oo, 

we have 

lim (max (1/X) ) = 0. 
The proof uses a standard branching process argument and is omitted. 



4 \i+2 



(!) 



The Bak-Sneppen model 



18 



4.1 Examples 

1. The original Bak-Sneppen model on the circle. Here Q = (P n )neN, 
where P n is the regular polygon with n vertices. Observe that for any k > n, 
we have M Pn (\2k/3]) = 0. For any k < n, we have 

m Pn ( [k/3\ ) > 2k/3, M Pn ( \2k/3] ) < 2. 
Then condition (|4.2|) holds for b > 60, and condition (|4.3J) holds for K = 3. 

2. The multi-dimensional Bak-Sneppen model. Let d > 2 be the 

dimension. We consider = (G n )neN, where for every n G N the set of 
vertices Vg„ = {1, ■ ■ ■ ,n} d , and with the usual nearest neighbour structure 
with periodic boundary conditions. Observe that for any k > n we have 
Mp n ( |~2A;/3] ) = 0. For any k < n, we have 

™p„(^/3J) ~ C!^, M Pn (\2k/3]) ~ c 2 k d ' 1 . 

Property holds for 6 > 30d, and condition (|Ojl holds for K = 2d+ 1. 

3. The Bak-Sneppen model on a tree. Choose d > 2 and consider 
Q = (Td(n)) n £N, where Td(n) is the regular d-ary tree with d offsprings 
at each vertex and n generations. Observe that for any k > 3n we have 
M Pn (\2k/3~\) = 0. For any k < 3n, we have 

m Pn ([k/3\) ~ c,d k /\ M Pn {\2k/3\) ~ c^ 2fc / 3 . 

One can check that property (|4.2j) holds for b > 30c?, and that condition (|4.3|) 
holds for K = d + 2. 

5 Proof of Theorem 14.11 

The proof of Theorem 14.11 essentially follows the proof of Theorem 11.11 For 
any finite G, one can associate a graphical representation GRq- The con- 
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struction of this graphical representation is essentially the same as before. 
Let {Hk}kGV G be a collection of independent homogeneous Poisson processes. 
For each process 11*. we perform the following procedure. At the j-th arrival 
r k)j of n fc , we draw a pair (£ G (k, r fcj ), rj G (k, r kjj ) ) , where £ G (k,T kJ ) is dis- 
tributed as the range set, fj G (k, r k j) as the duration of a typical r^.,- avalanche 
with origin at k. We draw arrows in V G X M + from (k, r k j) to (y, T kj j), for all 
y G icik, Tfcj). As before, we can define, for any A C V G , < t, s < oo, the 
processes ^'^(s) and ^'^(s), such that the monotonicity properties hold 

1o'"(si) <'i B '"(s 2 ), 0< Sl <s 2 . 
The range set ^g^, &) and the duration rj G (x, b) of a 6-avalanche with origin at 
x £ V G can be written as £ G (x,b) = ^^^'^(b) and r]N(x,b) = 7]^^^^^). 

The only place where we previously used the geometrical structure of G 
is Lemma 13.11 and the related definitions. We will give a new lemma for a 
collection of finite connected graphs. 

For any vertex x of G and any q > define Pc{q) as the probability that 
updating x and the neighbours in the configuration with all fitnesses above 
q results in a g-avalanche of range |Vg|- 

Lemma 5.1 Let G be a finite connected graph such that for some 60 < bo < 
oo, 

T (l) ™g( Ln/3J ) > log (m g { \2 ri /3] )) , uniformly imGN. 

Then there exists qooipo) G (0, 1), depending only on 6 such that for any 
q > <?oo0o) 

Pc{q) > 1/2, uniformly in x G G. 

PROOF: As in the proof of Lemma \i. II we will work with fitnesses defined on 
[0, oo), and we will update them according to the exponential distribution 
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with parameter 1, say. Recall that in the new setup a threshold b corresponds 
to the threshold q — 1 — e~ b in the old setup. 

We will have proved the lemma if we show that for any x G G, the process 
£g{ x , t) is explosive in the sense that there exists b^ G (0, oo), depending only 
on bo, such that for any i G N (essentially for i < i max = m&xp(y, z), because 

x,y 

for i > i max we have B G {ri) = V G ) 

P{BUn) C &(*,&«,)) > | + (|) m (5.2) 

Indeed, f)5.2j) implies that P^(6oo) > |. In order to achieve this, choose 
a constant &o satisfying the condition of the lemma. Define a converging 
sequence of thresholds bx, b 2 , b 3 , . . . as follows: 

h= h^ + dYbo, i>i, 

6oo = lim bi = 5b . 

Observe that due to the monotonicity property (|5.1jl . it suffices to prove that 
for all i G N, x G V G 

P(B G (ri) QtG(x,b i ))>l + (i) i+1 . (5-3) 
We proceed by induction. First note that for any x G V G 

P(flg(l)C £<?(*, &o)) = l. (5.4) 
Next, suppose that (|5.3jl holds for some i G N and all a; G Vq. Observe that 



< Vz G 5g(r2r</3l), 3y G 5^(^/3]), 3r w - G II, n Mm) 
such that B y G (ri) C £ G (y,T yJ ) 

implies 

{ Bg(r i+1 ) Cfe(x,6 i+1 ) }• 



} (5-5) 
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Indeed, if a G B G (r i+ i) \ B G (ri) then there exists z = z(a) G S x (\2ri/3\) 
such that p(a,z) + p(z,x) < r i+1 , and hence p(a,z) < |~2r.j/3~|. Then if 
there exists r G |J II,, D [bi,b i+ i) such that B v G {ri) C ^g(z/) t ) then 

p( a ,y) < p( a iZ) + p(z,y) < ri, and hence a G £g(?/,t). See Figure |U for an 
illustration of (|5.5jl . Hence to finish the inductive step, it suffices to show 




Figure 4: Illustration of (|5.5jl . 



that, uniformly in i G N, 



P 



3z G ££([2^/3]) such that 
Vy G Sg(Lri/3j), Vr^ eU y n M*+i), 

B y G { ri )%Uy,^) 



\ 



<(§) 



l\i+2 



(5.6) 



Since the events B v G {ri) C ^(y, independent and since (due to the 

monotonicity property) for any y G Vq and r y j > 



P KN c eo(y,r WJ )) < P (B G { n ) c &(y,6i)J < I 



2' 
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for any z G ( [2r j /3~| ) the points 

U \r y ,j e Ii y ; B y G (n) C | G (|/, r tf j) } (5.7) 

constitute a thinning of the Poisson process 1J U y with deleting prob- 

yeB5(Lri/3j) 

ability at most 1/2. Thus the points in ()5.7j) contain a Poisson process of 
intensity at least 

- ■ (intensity of (J n y j > - ■ m G (Lr;/3j). 

yes&(Ln/3j) 

Observe that the event in ()5.6|) implies that for some z G SQ(\2ri/3~\), the 
process (15. 7j) has no arrivals between time 5j and &i+i, a time interval of 
length bo (|)\ The last event has probability at most 

exp|-|m G (Lri/3j)6 (f)*} , 

uniformly in z G Sq ( |~2r j /3] ) . Thus we can estimate the probability in (j5.6J) 
by 

M G ([2r l /3l)exp{-|m G (Lr i /3j)(|) J }. 
Split the above expression in two terms 

(M G ([2r,/3l)exp{-^m G (Lr l /3j)(|) J })exp{-|m G (Lr i /3j)(|) i }. 

The first term is less than or equal to 1, under the conditions of the lemma. 
For the second term, we write 

exp {- b -fm G ([r,/3|) H) 1 } < exp {-| (1 + Lr./Sj) (§)*} 
< exp {~ b i (§ (|) 4 } < exp (If) 4 } < (I) 4+2 , i G N, 

since 60 > 60. So we have (|5.fi|) and the proof is complete. 

□ 
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